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Aphaenogaster longiceps, the funnel ant, is widespread in southern Australia, The genus is 
generally regarded as being more common on, or conlined tv, sandy soils ranging from the 
inland to the coastal forests. In this paper we model the effect of site, soil and vegetation 
characteristics on the probability of occurrence of nests of A. longiceps. The relau ve strength 
of the environmental variables in predicting the presence of nests of A. longiceps within a 
20ha site in SE New South Wales is described. Aspect and topographie position are belter 
predictors of nest occurrence than are surface soit texture (pereen! gravel and percent sand) 
or vegetative characteristics. The models are used to predict the likely response of A. 
longiceps to foresi habnat fragmentanon in eucalypl furest in SE New South Wales. 
(Hymenoptera, Formicidae, funnel ant, Aphaenogaster longiceps, spéecles distribution, 
statistical modelling, logistic regression, habitat fragmentation, Wog Wog, Australia. 
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The impuct of forest habitat fragmentation on 
the ground dwelling arthropod fauna of a 
cucalypt forest is currently being studied at Wog 
Wog, SE New South Wales (Margules, 1992). 
The primary objectives of this experiment are to 
test predictions that flow from MacArthur & 
Wilson's (1967) equilibrium theory of island 
biogeography. The theoretical background, cx- 
perimental design, sampling stratification and 
theoretical predictions are described by Margules 
(1992), In brief, there are two predictions of im- 
portance: the first is that habitat fragmentation 
reduces diversity and second that the reduction in 
diversity is dependent upon the fragment size 
(Margules, 1992). The major emphasis of the 
experiment is on the ground dwelling arthropod 
fauna. Herbage and understorey layers are also 
being monitored (Austin & Nicholls, 1988; Mar- 
gules, 1992) to document their dynamics follow- 
ing habitat fragmentation and for the polential 
influence they might have on the distribution and 
abundance of the sampled fauna. 

Funnel ants, Aphaenogasterspp., are known tor 
their habit of creating widespread surface soil 
disturbance in association with their nest entran- 
ces (Sloane & Sloane, 1964; Saunders, 1967; 
Andersen, 1991). They can be a conspicuous 
component of the ground dwelling fauna and can 
influence the dynamies of the ground layer 
vegetation (Saunders, 1967). The presence of the 


nests of A. longiceps was recorded as part of the 
characterisation of the permanent monitoring 
sites in the Wog Wog experiment. The nests did 
not seem to he randomly distributed with respect 
10 other site environmental variables. 

The preference of Aphaenogaster spp. for 
selected habitats has been noted belore 
(Saunders, 1967; Andersen, 1991); itis generally 
observed that Aphaenogaster spp. prefer sandy 
soils, Although there has been only limited work 
on habitat preferences of Australian ants there is 
clear evidence that ants can show affinities for 
selected soil types (Greenslade, 1976, 1987). 
aspects (Greenslade, 1985) and can partition the 
habitat in terms of rainfall gradients (Greenslade, 
1974, 1976, 1987). 

The objective of the present paper is todevelop 
and describe correlative models relating the 
presence of nests of A. longiceps to a range of 
environmental variables that potentially repre- 
sent different sources of contro) of the distribu- 
tion of the species. These models will be used to 
predict the likely impact of forest fragmentation 
on this conspicuous component of the ground 
dwelling arthropod fauna. 


METHODS 


In 198? an environmental survey of the 144 
permanent sampling pomis within the Wog Wog 


TABLE 1. Variables recorded for each permanent 
monitoring site grouped into three broad categories. 
Categorical variables or factors are indicated by the 
number of levels within brackets. Soil colour was 
recorded using the Munsell Color Company (1971) 
hue, value, chroma polation; electrical conductivity 
(EC) in a 1:5 water solution. Organic carbon is ab- 


breviated as OC_ 

Topographi 
Slope |Ahorizon |% Bare ground 
[Aspect | Sand Litter cover | 
[Water Erosion(2) |__| Eucalyptus, other) 
|  [Bhorizon __|% Grass cover_| 


% Herb cover 
cover 


moss, lichens 
% Cover fems 
% Canopy cover 


(Eucalyptus, 


Acacia) 


| Organic Carbon | O 


habitat fragmentation experiment was carried out 
just prior to the clearing of the surrounding forest. 
Variables were recorded in three broud classes: 
topographic, edaphic and vegetational (Table 1). 
Each permanent monitoring site consists of two 
3X3m quadrats. separated by distances that 
varied from two to about four metres, on which 
most vegetational yariables were recorded. The 
only vegetational characteristics not measured on 
these quadrats Were the canopy cover of Eucalyp- 
Ius spp. and Acacia spp., which were estimated 
for a 10% 10m quadrat centred on a soil auger 
hole placed between the two 33m quadrats. 
Edaphic data were collected from this auger hole 
and associated samples submitted for laboratory 
analysis. Topographic data were recorded from 
the same location as the edaphic data. Presence 
or absence of nests of A. longiceps were recorded 
for each 3 X 3m quadrat. 

The relationship between the presence or ab- 
sence of nests of A. longiceps and the different 
suite of environmental variables was determined 
by fitting a statistical model — a multivariate 
logistic regression, from the class of regression 
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models known as generalised linear models (Mc- 
Cullagh & Nelder, 1989). Numerous examples 
of the application of logistic regression models 
to demonstrate the relationship between the 
presence or absence of a species and a suite of 
environmental variables have been published for 
both plant species (Austin et al., 1983, 1984, 
1990: Nicholls, 1989) and animal species (Adler 
& Wilson, 1985; Lindenmayer et al., 199 1a, 
1991b). Models were developed using the for- 
ward stepwise variable selection strategy out- 
lined by Nicholls (1989) for predicting the 
probability of occurrence of Eucalyptus radiata 
as a function of aset of climatic, lopographic and 
geological variables. The presence of a 
monotonic response as opposed to a unimodal 
response in each continuous variable was tested 
by fitting and comparing both first and second 
order polynomial functions. 

A separate model was developed to estimate 
the probability of nest occurrence as a function 
of the variables belonging to cach of the classes 
listed in Table 1. This permits an intuitive ap- 
proach to evaluating the relative strength of the 
three models while permitting one to look at the 
variables that are correlated with the presence of 
Aphaenogaster nests. A more formal approach 
to this relative evaluation would be to use the 
forward selection procedure across all variables 
from all three classes. The usual set of regression 
diagnostics, techniques developed to assist with 
the evaluation of the assumptions. implicit in 
fitting regression models to data, were applied ta 
the final models to assess observations not well 
fitted by the model (outliers) and observations 
with undue or potential influence on the 
parameter estimation (Nicholls, 1989; Hosmer 
& Lemeshow, 1989; Collett, 1991) 


MODELS FOR INDIVIDUAL CLASSES 
OF VARIABLES 


TOPOGRAPHIC VARIABLES 

A number of variables showed significant 
relationships when tested singly. Of these the 
following were the most important: aspect, the 
expenmental design factor habitat (slope or 
drainage line), catchment size and landform 
(Table 2). The variables — habitat, landform and 
catchment size — produced a significant im- 
provement when added to a model with aspect 
(Table 3). The variables slope and landform 
produced significant improvement in a model 
including aspect and habitat. Of these, slope was 
the most important and was added to the model. 
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TABLE 2. Changes in deviance (a measure of the lack of fit used to assess logistic regression models fitted to 
binary data) when the topographic variables were fitted singly to the null model. This table represents the results 
of the first pass of the forward stepwise fitting procedure. 


Null 287 213.09 


Subsequent testing of the remaining vanables 
produced no further significant reduction in the 
deviance. Thus a three variable model proved to 
be an adequate model. The regression coeffi- 
cients for this mode] are given in Table 4 and the 
model displayed in Fig. 1. Interpretation of the 
regression coefficients is as follows: the positive 
cocfficicnt for aspect indicates an increasing 
probability of finding a nest as the aspect in- 
creases towards 360°, the negative coefficient for 
slope reflects a declining probability of nest oc- 
currence on flatter sites compared to siles on 
slopes. The habitat coefficient is the difference 
between the two predicted surfaces defined by 
aspect and slope and reflects a 3.5 times decline 
in the ratio of the odds (presence/absence) of a 
prediction for a drainage site compared with a 


20 
15 
a 
& 10 
D 
5 
0 
100 150 200 250 300 350 
Aspect 


afo | Devine | af. | Deviance | =P 


i eet Bs ee ee 
| 1.097 


[Siope | 6e | 2.68 
Slope+Slope* 285 211.59 
Aspect 286 181.80 
sas i 
279 180.51 
286 194.55 
254 193.19 
210.02 


31,294 
I 3.353 0.067 

32,582 0.001 
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comparable slope site. Figure 1 shows the dis- 
tribution of the observations in the multidimen- 
sional space defined by slope and aspect. The 
lines show wherc in that two dimensional space 
the mode! predicts selected probabilities of nest 
occurrence for given combinations of slope and 
aspect. Note that a predicted probability of 0.1 
implies that for ten sites (for example slopes) 
with the same slope and aspect, one would ex- 
pect one site to have Aphaenogaster nests; in 
contrast for a predicted probability of 0.5, five of 
the ten sites would be expected to have nests. 


EDAPHIC VARIABLES 

Textural variables, percent sand, gravelor elay 
for the Al horizon or the B horizon showed a 
significant relationship with the presence of A. 
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FIG. 1. Observed and predicted probabilities of occurrence of nests of A. longiceps as a function of aspect, 
expressed as degrees from north, for sues classified in the ficid as either slopes (a) or drainage lines (b). Sites 
are shown as solid symbols where nests were present and as open symbols where absent. Continuous lines 
show the predicted probability of nest occurrence, 
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FIG 2, Distribution sites in the environmental space 
defined by percent gravel and percent sand. Sites are 
shown as solid symbols where nests were present and 
as open symbols where absent. Continuous lines 
show the predicted probability of nest occurrence. 


longiceps nests (Table 5). The colour, expressed 
as the chroma (Munsell Color Company, 1971), 
of the A2 and B horizons showed less pronounced 
but significant relationships (Table 5). The final 
model was a quadratic funetion of the percent 
gravel of the A} horizon plus a linear function of 
the percent sand of the A] horizon (Table 6). For 
any given values of percenl sand the maximum 

redicted probability of occurrence will be at 
about 15-16% gravel content, The predicted sur- 
face and observed data are shown as a function of 
these two variables (Fig. 2). 


VEGETATIONAL VARIABLES 

The vegetational characteristics fal] naturally 
into three groups: litter, understorey components 
and the canopy. Significant relationships were 
found for litter and understorey components but 
not for canopy components (Table 7). Percent 
bare ground was the most significant litter com- 
ponent whereas percent grass cover was the most 
important single variable of the understorey 
group (Table 7). The selected model was a quad- 
ratic function of percent bare ground, linear func- 
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tion of percent moss cover and a linear function 
of percent grass cover (Table 8). Maximum 
probabilities are predicted for percent bare 
ground around 18%; the probability of occur- 
rence declines with increasing cover ot grass and 
moss. For moss cover the rate of decline of the 
odds ratio is about 1.6 times that for grass cover. 


COMPARISON OF THE THREE MODELS 


The three models can be compared informally 
in terms of the overall measure of laek of fit. For 
generalized linear models, this measure is usual- 
ly referred to as the deviunce (McCullagh & 
Nelder, 1989). Although for many generalized 
linear models the behaviour of the devianee is 
understood for models based on binary (that is 
presence or absence) data (such as used here) the 
behaviour is known not to have expected 
asymptotic properties. For this reason the use of 
the deviance as a measure of model adequacy is 
not recommended (McCullagh & Nelder, 1989). 
However, we can use it to rank the three final 
models; the topographic model had a residual 
deviance of 152 with 284 degree of freedom; the 
vegetation model a deviance of 168 with 283 
degrces of freedom; and the soil-based model a 
deviance of 190 with 284 degrecs of freedom. 
The topographic-based model offers the greatest 
explanation of the observed variation in distribu- 
tion of nests of A. longiceps and the vegetation 
model the least explanation, 


DISCUSSION 


Despite the rarity of occurrence of nests — 
present on 35 of the 288 quadrats (which were 
located on 25 of the 144 sites) — a wide range 
of cnvironmental variables drawn from the three 
classcs have significant relationships with the 
occurrence of A. longiceps nests (Tables 2, 5 & 
7). This result reflects the likely collineanty of 
many of the environmental variables within the 
sample of sites used for monitoring the impact of 
fragmentation, Despite this, itis possible to draw 
some conclusions about the influence of en- 
vironment on the distribution of A. longiceps 
within the local region of the Wog Wog habitat 
fragmentation experiment. 

The site or topographic-based model 
demonstrates a clear preference of species for 
warm exposed (north-west facing) slopes (Fig. 
1), In addition, the preference for sandy soils on 
warm slopes suggests that this species avoids wet 
conditions. As noted by Margules (1992) the 
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TABLE 3. Changes in deviance (a measure of thc lack of fit used to assess logistic regression models fitted to 
binary data) when the topographic variables were added singly to the model containing the linear function on 
aspect. This table represents the results of thc second pass of the forward stepwise fitting procedure. 


| Variable | df | Deviance | df | Devine | P | 
aseet f so o mM 
ae a a ee 
reer | zea [ana apey so 
Suet ae ae le ee | ee 


[+Landform | 278 | a | e o aa |<. 
Hasia [as ooo o aon o01 
l+Catchmentsize | 283 | 16036 | 3 | ts | <o | 
Fenian [aes | neas f o asa oos 
[+ Water Erosion | 285 | 1699 |  ı | ass | 0o07 | 


TABLE 4. Regression coefficients and approximate standard crrors for the final model relating probability of 
occurrence of nests of A. longiceps with selected topographic and site variables. 


| Variable |Parameterestimate| Standarderror | tvalue | P | 
1.598 
i 


0.03468 


-0.2407 0.0839 2.87 0.005 


TABLE 5. Changes in deviance (a measure of the lack of fit used to assess logistic regression models fitted to 
binary data) when the edaphic variables werc fitted singly to the null model, Textural components are 
abbreviated, sand = S, clay = C, and gravel = G. All continuous variables were fitted as both first and second 
order polynomial functions. Where no second order function is shown below it may be assumed that the change 
in deviance due to the addition of the quadratic term was not significant and was less than 2. 
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TABLE 6. Regression coefficients and approximate standard errors for the final model relating probability of 
occurrence of nests of A. longiceps with selected edaphic variables. Residual deviance for this model was 
190.05 with 284 degrees of freedom. Edaphie variables are abbreviated, % gravel as %G and % sand as %S. 


Standarderror | tvaue č | O P | 
Intercept -5.9428 1.8491 


Al %G 0.2639 0.0620 <0,001 
Al %G + pG? -0.00844 0.00305 297 0.007 
Al %S 0.0496 0.0258 


TABLE 7, Changes in deviance (a measure of the lack of fit used to assess logistic regression models fitted to 
binary data) when the vegetation variables were fitted singly to the null model. All continuous vanables were 
fitted as both first and second order polynomial functions. Where no second order polynomial function is shown 
below it may be assumed that the change in deviance duc to the addition of the quadratic terin was not significant 
and was less than 2. 


Null 213.10 


% Bare grnd 286 
Gh Bare emd” 285 188.94 


[Variable | af | Deviance | af | Deviance | P | 


I 12.49 <0.001 


TABLE 8. Regression coefficients and approximate standard errors for the final model rclating probability of 
occurrence of nests of A. longiceps with selected vegetation variables. 


-1.2785 
0.3249 
-0.0110 
0.0972 


intercept 


% Bare ground? 
% Moss 
% grass 


experimental site was selected such that the 
layout of the remnants was on a predominantly 
south western facing slope. Despite this, the local 
relief within the experimental site causes aspects 
for the 144 monitoring sites to range from 107° to 
325° with the majority between [85° and 240°. 
There is little difference between those monitor- 
ing points classified as slopes and those as 
drainage lines in terms of the range of aspects nor 
in terms of their slope. There is little correlation 
between site aspect and slope (Fig. 1). The impact 
of fragmentation of the forest will not change the 
aspect or slope of the monitoring sites although 
the exposure of sites can increase due to the loss 
of surrounding forest canopy This can he most 


0.4434 
0.1254 
0.0072 
0.0428 


-0.0596 0.0267 


0,010 
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0.026 


noticeable for sites close to the edges of the 
remnants and may change the range of effective 
exposure, increasing it for westerly aspects and 
perhaps decreasing it for southerly aspects. 
Given the preference of A. longiceps for nesting 
on westerly to north-westerly (=warm) sites and 
the potential for increased exposure following 
forest canopy fragmentation, the response of this 
ant to the direct impact of fragmentation will be 
to maintain or perhaps to increase its current rate 
of site occupancy. 

The strong influence of surface texture on the 
probability of nest occurrence is consistent with 
the usual statements about the preference of 
Aphaenogaster spp. for sandy soils (Sloane & 
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Sloane, 1964; Saunders, 1967; Andersen, 1991). 
Although the best model includes the Al horizon 
sand content, it also contains the gravel content 
of that layer, a point not noted in the literature. 
Neither of these two variables are likely to change 
dramatically due to the fragmentation treatment- 
The question of collineanty within the data 
matrix is important because the soils with high 
clay content tend to dominate the gullies. Total 
separation of the influence of surface texture from 
that of topography might not be possible without 
more extensive sampling. Although a formal test 
of the difference between the two mode!}s has not 
been undertaken, the site-based model appears to 
offer better predictive ability than the soil-based 
model. 

The interpretation of the final vegetation-based 
model is less clear because more variables are 
included and these show more correlation than 
one would bike when undertaking this type of 
modelling exercise because of the difficulty of 
separating the influence of individual variables 
when there is high collineanty. Also, these van- 
ables can themselves respond to the impact of 
fragmentation of the forest. Increased exposure 
— particularly around the edges and perhaps over 
all of the small remnants — may result in a loss 
of moss cover and a reduction in grass cover. This 
could lead to an increase in bare ground, one of 
the important predictor variables in this model. 
Such changes, if they occur, may lead to an 
increase in favourable nest sites because the prob- 
ability of occurrence is negatively correlated with 
both moss and grass caver. 

It should be noted that there are substantial 
differences in the distribution of moss cover on 
the slopes in contrast to the drainage lines. The 
inclusion of moss cover in the model as a predic- 
tor might be the result of this difference, Frequent 
and high moss cover is characteristic of drainage 
lines; as such, it may be acting as a surrogate for 
the habitat variable noted as an important predic- 
tor in the site-based model rather than the anis 
responding directly to the extent of moss cover. 

On the basis of the three independently fitted 
models, the dominant influences on the distribu- 
tion of nests of A. longiceps appear to be those 
that are not likely to change following the clear- 
ing of the surrounding forest to create habitat 
remnants. This suggests that the fragmentation 
treatment imposed on this eucalypt forest ts un- 
likely to have a direct impact on the distribution 
of A longiceps nests. This is not im conflict wath 
the concern expressed by Sloane & Sloane (1964) 
with the impact of habitat loss on A. longiceps. 
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The prediction that fragmentation will not 
result in a reduced rate of site occupancy will be 
tested in two ways and the results reported else- 
where. A re-survey of the quadrats is planned for 
the spring or summer of 1993/94, This would be 
directly comparable to the 1987 survey and nest 
loss could be determined and related to the isola- 
tion treatment, habitat remnant size and position 
within the remnant, A second test of the predic- 
tion will be based on changes in distribution and 
relative population size of A. longiceps caught in 
pitfall traps operated quarterly for the duration 
of the experiment (Margules, 1992). 
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